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ABSTRACT 


The wind motions responsible for the shearing of sodium 
vapor trails ejected from rockets in the 70 to 140 km region 
of the upper atmosphere are subjected to an analysis based on 

generally accepted theories of hydrodynamic turbulence. The 

region from 80 to 100 km is of particular interest in that here 

the predictions of shear turbulence theory are well substantiated. 

The energy spectrum of the height shear is found to follow the 4/3 

power law proposed by Tchen, and is associated with a vertical 

correlation distance of approximately 6 km. The existence of an 

isotropic inertial region of maximum scale 3 km, previously 

indicated by analysis of meteor data, is confirmed. The vertical 

scale of the turbulent eddies is found to be the atmospheric 

pressure scale height, a phenomenon which has been observed by 

others, but which, as yet, has no satisfactory explanation. 



1 . Turbulence Theory 


The complete development of the relationships used in the 
analysis of wind shear is beyond the scope of this work; the 
following, used in conjunction with the references quoted, 
should provide an adequate background for consideration of the 
subsequent analysis. 

1 . 1 Energy Spectrum Analysis 

If there exists in a turbulent flow field a range of 
scales which receive energy from larger scale motions and pass 
it on undiminished to smaller scale motions, then, for this so- 
called inertial (nondissipative) range of scales, the only 
form of energy spectrum function dimensionally possible is 

E(k) = ae 2/3 k " 5 ^ 3 (Kolmogorof f , 1941) 

where e is the rate at which the turbulent energy is received by 
(and leaves) the inertial range of scales, k is the wavenumber 
vector corresponding to the real space scale r, and a is an 
absolute constant of order unity, 

Batchelor (1953) has shown that, for such an inertial 
region which also possesses the property of isotropy, the fluid 
velocity differences measured as a function of the separation r 


follow the relation 
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[u(x)-u(x + r) = 4.82a(er) . ....... .1 

In real space, the energy spectrum function E(r) is 
defined by 

2 

E(r) = [u(x) - u(x + r) ] o .. .2 

Tchen (1954) has considered an otherwise isotropic 
region subjected to a mean shear, and finds that equation 1 is 
modified, becoming 


[u(x) - u(x + r) ] = 


OOOCCOOOOO*QO< 


where a involves a, e, and the mean gradients Thus, for what 
may be termed shear turbulence, in real space 


E(r) 


o*«*o«««oooooooo*c 


1 . 2 Correlation Analysis 

An energy spectrum function equivalent to that based on 
velocity differences can be formulated from the lateral or 
longitudinal velocity correlations defined as 


f (r) 


u f (x) u f (x + r) 
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and 


u (x) u (x + r) 

g(r) = ^ 


oo**o*ooo»**ooo< 


u 


n 


where u f (x), u f (x + r) , u^Cx), u n (x + r) are the turbulent 

components of the velocity at two points x and x + r respectively, 

measured parallel (suffix f) and normal (suffix n) to the vector 

2 2 2 

separation r„ In isotropic turbulence, u f “ u n = U Q > the 
velocity characteristic of the energy bearing eddies of the 
Kolmogoroff spectrum,. 

Introduction to the equation of continuity for 
incompressible fluids leads to the relation 


g(r) f (r) i - 1/2 r ^ ~ ooo*»o»o«ooooooooo6 

(von Karman and Howarth, 1938) 

between the functions f and g, or, if the turbulence is isotropic 
in two dimensions only, 


g(r) 


f (r) 


+ 



ooooavco* 


,7 


These functions, f and g, may be called Eulerian space correlation 
functions, and either may be denoted by R(r) „ If the range of 
scales under observation is inertial, then R(r) must depend only on 
e, and the only form dimensionally possible is 
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u Q [1 - R(r) ] 


2/3 2/3 

e r 


i .e „ 1 - R(r) ~ r 

(for any given flow field, u Q and e are constant). 
Introduction of the relations 6 and 7 above gives, 


with c a constant 


f = 1 - cr' 


. . .8 


g = 1 - 4/3 cr /J t . . . . . .9 
for three dimensional isotropy 

-o _ t c /o „ in 


f = 1 - 5/3 cr 


...10 


for isotropy in two dimensions only. 

Since the variations of the correlation difference 
function [l-R(r) ] with r is the same as that of E(r) in equation 1, 
we may suppose that the dependence of [l-R(r) ] on r will be 

similarly modified in the presence of a mean shear 
j « n r> s i -4/3 , , 


[l-R(r) ] ~ r 


OOVVOOOOO. 


If, in fact, equation 11 is pertinent, then the 

relationships between f, g, and r (equations 8, 9, 10) will now become 

4/3 

1 = 1 — ar 

and g = 1 - 5/3 ar^^ 

or g = 1 - 7/3 ar 4/3 

An indication of the degree of isotropy can be obtained 


by considering the ratio 


1 - f 
T-=”g 
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For two dimensional isotropy 

S = 0.60 without mean shear 
and S = 0,43 with mean shear 

For three dimensional isotropy 

S = 0,75 without mean shear 
and S = 0,60 with mean shear 

The relative importance of mean shear will be indicated by 
the form of either the energy spectrum function E(r) , or the 
correlation difference function [1 - R(r)]„ 

2 . The Practical Application of Turbulence Theory 
2.1 The Mean Wind Profile 

In applying the relations developed above to the wind 
vectors measured by means of sodium vapor trails ejected as a 
tracer into the upper atmosphere, the relative importance of the 
mean motion must not be overlooked. The velocities used in the 
energy spectrum and correlation analyses of the previous section 
must be the turbulent velocities, or departures from the mean 
motion. In normal correlation analysis, the mean of a set of 
observations is usually taken as the mean value of the measured 
data, which, when applied to wind height shear data, would be 
tantamount to the assumption of a mean wind profile which is 
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constant with height. Such a profile is the exception rather 
than the rule in meteorological phenomena. 

Practically any attempt to prescribe a mean wind profile 
to the measured winds will be subjective to a certain extent. 

If one uses a polynomial fdt, for example, it must be truncated 
before it can assume any of the features of the measured profile 
which are due to the turbulent motions present. From experience 
based on the measurement of winds in the 75 to 105 km region by 
means of radio reflections from meteor trails (Elford, 1958, 1964), 
a quadratic change with height of the mean wind over any given 20 km 
interval should best describe the contribution of the mean motion 
without destroying any of the characteristics of the turbulent 
flow field. In the present analysis, the windspeed/azimuth data 
are converted to zonal and meridional components, and a polynomial 
profile of order Z + 1, where 

„ _ Total height range covered by data(km) 

z - • — — j 

20 

is fitted to each. The profiles thus determined are subtracted 
from the relevant measured profiles to give zonal and meridional 
turbulent velocities. 
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2 „2 The Determination of E ( Ah) 

The available sodium trail data lists wind speed and 
azimuth against height, and usually involves irregularly spaced 
height intervals. The sampling irregularity exists for two 
reasons : 

a) there is an occasional difficulty in absolutely 
identifying the same point on the trail in consecutive photographs; 

b) the wind profile between consecutive observational 
heights is linear. This is usually obvious from the photographs, 
and can easily be allowed for in subsequent analysis. 

Whereas equal height interval sampling is not absolutely 
necessary for subsequent reduction, it does simplify the analysis, 
and so linear interpolation between the listed data points is used 
to provide a profile with data points spaced 0,2 km in height. 

The energy spectrum function is computed as 

1 N 2 

E ( Ah^) ^ 2 (u^i-u^) .oo.oooo 12 

k= 1 

in which 

N = SCHg-Hj) - i 

where i = 1, 2, 3....N 
such that Ah^ = 0.2 km. 

Ah ^ = 0„4 km 


etc o 
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and (in km) are the lower and upper bounds respectively 

of the region for which E(Ah) is being determined. Such 
partitioning of the height range is necessary since there is 
considerable variation in the characteristics of the flow over 
the total height range sampled (usually some 70 to 200 km) . 

Energy spectrum functions may be calculated using 

a) zonal turbulent velocities 

b) meridional turbulent velocities 

and c) turbulent windspeed. 

The turbulent windspeed w is defined here as 

w = measured windspeed - ((mean zonal wind/ + 

2 1/2 

(mean meridional wind) ) ' 

2 0 4 Correlation Analysis 

The two correlation functions f and g defined in 1.2 
refer to turbulent velocity components measured parallel and 
normal to the separation vector. Since vertical velocities in this 
region of the upper atmosphere are so much less than the associated 
horizontal components, the magnitude of the vertical component 
cannot be determined from sodium trail photographs. (Most workers 
in this field consider the upper limit for mean plus random vertical 
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motions to be some 10 metres/sec) . However, we may redefine 

f and g in terms of the orthogonal zonal and meridional flow 

fields, with a view to investigating possible isotropy. This has 

been done by meteorologists in the past, with at least partial 

success (see, for example, Hutchings (1955)). 

2 2 

The normalizing factors u^ and u Q are best estimated 
by the standard correlation function definition. The correlation 
functions f and g then becomes 


f (Ah ± ) 


N 

£ 

k=l 


u 


k+i 


u 


k 
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where u are the zonal turbulent wind velocities and 


g( Ahu) 


L v * 


[i 

k=l 



N 

I 

k=l 



~T7Z 

] 


J 


15 


where v are meridional turbulent wind velocities, and N and i are 
as previously defined for equation 12. 
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3 . Preliminary Results 

The analysis of Section 2 has been applied to data obtained 
from a sodium trail release over the Eglin Air Force Base, 

Florida (29.6°N, 86.6°W) at 1910 CST on May 21, 1963 (Edwards 
et al, 1963). In this experiment, wind speed and azimuth were 
obtained over a height range of 69 to 140 km. 

3.1 The Mean Wind Profile 

Since the data covers the height range from 69 to 140 km, 
polynomials of order 4 are fitted to the zonal and meridional 
measured profiles. These yield mean zonal and meridional profiles 

u = 38.6 - 126h - 244h 2 + 174h 3 + 205h 4 
mean 

v = -30.4 + 19. 8h + 148h 2 - 16. 8h 3 - 125h 4 
me an 

where u, v are in metres/sec and h is the normalized height 
given by 

h = (2z - z , - z )(z. -z ) 

min max mm max 
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where z is the height variable 

z the maximum and 
max 

z min the minimum heights of the available data, 
all heights being in kilometers „ 

Normalization of the height range stabilizes the 
least squares fitting process, and makes the relative importance 
of the individual terms of the fitted polynomials more obvious 
than is the case when the mean velocities are expressed as power 
series in the actual height z. Results are plotted in Figs.l and 2 „ 
The meteorological significance of these profiles, in 
particular the reversal of both the zonal and meridional components 
above 110 km, cannot be evaluated from consideration of this 
single firing,, 

3 .2 The Turbulent Wind Profile 

As mentioned in Section 2, the 70 to 200 km height range 
covers a number of characteristically different regions. The wind 
motions observed below approximately 105 km indicate the presence 
of small-scale structure, while those above 110 km do ..not appear 
to be at all turbulent, even though vertical shear is present. 

In the results presented here, discussion is confined to the 
consideration of the region from 80 to 100 km, which has been found 
to be representative of a turbulent region which can be adequately 
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described by available statistical theories of hydrodynamic 
turbulence (Blamont and Jager, 1961; Zimmerman, 1962; Roper, 1962) „ 
The measured zonal and meridional profiles, and the 
deviations from the mean wind for the height range 80 to 100 km 
are plotted in Figures 3 to 6. Whereas the immediately obvious 
wavelike nature of the turbulent profile would suggest a wave 
theory approach as likely to be the most profitable for consideration 
of the wind motions in this region, applications of available wave 
theories (e.g., Hines, 1959) have not produced consistent results „ 

The possibility of generation of turbulence in the 80 to 100 km 
region by vertically propagating gravity waves has been proposed 
by Hines (1963) . The purpose of the present work, however, is not 
to determine the source of the turbulent energy, but rather to 
substantiate the evidence that the observed shears characterize 
a region of hydrodynamic turbulence „ 

3 o 3 The Energy Spectrum Functions 

The energy spectrum functions for the zonal and meridional 
turbulent wind profiles for the 80 to 100 km region have been 
computed using the methods described in Sections 2.2 and 2.3. 

As can be seen from Figures 7 and 8, the velocity difference 
spectrum E(Ah) and the correlation difference function [l-f(Ah)] 
for the zonal turbulent velocities are completely equivalents 
This is not surprising, since both functions are solutions of the 
same equation 
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F = a Ah m 

(see Sections Id and 1„2). 

Similarly for the meridional spectra of Figures 9 and 10 „ 

The slope m ofthe log log plots of both the zonal and meridional 
E(Ah) (or [l-R(Ah) ]) spectrum functions against height difference 
Ah is constant at 4/3 for separations Ah up to 3 km, indicative 
of an isotropic region subject to a mean wind shear (Tchen, 1954) „ 

3 „4 Isotropy 

Since the energy spectrum functions and the correlation 
difference functions all follow an established (shear) law at small 
scales in the height range from 80 to 100 km, it is possible to 
determine the nature of the isotropy at these scales in this- 
region . 

Figure 13, curve A, is a plot of the variation of 

s = i=£ 

i-g 

against Ah« 

In the region up to a scale of 3„5 km, S has a value of 
approximately 0„7, indicating an isotropy lying somewhere between 
two and three dimensional. The increase in S at scales greater 
than approximately 3 km is due to the breakdown at this scale of 
the 4/3 power law in both zonal and meridional spectrum 


functions . 
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The function S has been found to be quite sensitive 
to variation in the parameters specifying the mean wind profiles. 

For example, if, instead of the fourth order polynomial fits to 
the total height range of the zonal and meridional data, linear 
mean wind profiles are fitted over the height range 80 to 100 
km only, subsequent spectrum analysis yields the S function of 
curve B of Figure 13. This curve would indicate the existence 
of three dimensional isotropy for scales up to 3 km. It is 
possible that the 4th order polynomial fits are, in fact, attributing 
a small fraction of the random wind variations to the mean motion. 
However, in order to correlate spectra obtained for different height 
strata (as is done in the next section) , a continuous mean profile 
is required, and has therefore been used throughout the analysis. 

3.5 V ertical Scale 

The vertical scale associated with the turbulent wind 
structure is conveniently defined by the Ah corresponding to the 
maximum value of either E(Ah) or [l-R(Ah)]. As can be seen from 
Figures 7, 9 and 11, which are plotted using zonal, meridional 
and windspeed turbulent velocities respectively, the vertical 
correlation distance as defined above is not the same for each 
component, i.e., the wind motions are not isotropic at the maximum 
of E ( Ah) . However, for the purpose of comparison with stratospheric and 
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lower mesospheric data, where meridional winds are, for the 
most part, negligible, consideration need only be given to the 
zonal energy spectrum. 

The variation of vertical scale with height in the 
stratosphere has been determined by Webb (1964) from a series of 
Robin soundings at Eglin. His results indicate an exponential 
increase of the vertical scale with height, from approximately 
800 metres at 35 km to 2 km at 55 km. He has extrapolated this 
exponential to a height of 90 km, and finds a vertical scale of 
6 km which is the value determined for this height by Greenhow (1959) 
from a correlation analysis of wind shears determined by means of 
radio reflections from meteor trails at Jodrell Bank (53°N) . 

In order to investigate the change with height of the 
vertical correlation distance in the 80 to 130 km region, the 
sodium trail data was subjected to a stepwide analysis, starting 
with the 75 to 95 km height range, and proceeding via the 80 to 
100, 85 to 105, etc. ranges to 115 to 135 km. The maxima of the 
resultant 'E(Ah) curves were then plotted as the vertical scales 
at the midpoints of the respective height ranges. The results 
are shown in Figure 14, together with Greenhow 1 s determination, 
and the 7.8 km at 94 km (October, 1961) determined from radio 
meteor trail shears at Adelaids (35°S) (Roper, 1962) . The dashed 
line is Webb’s extrapolated variation. While there is excellent 
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agreement with Webb’s proposed exponential increase, there 
is also good agreement between the vertical scales measured 
at various altitudes for the Eglin firing, and the atmospheric 
pressure scale height (also shown in Figure 14) „ This phenomenon 
has also been observed by others (e,g., Zimmerman (1964) )„ As 
yet, no satisfactory explanation for such a dependence has been 
proposed . 

Above 125 km, the magnitudes of the deviations of the 
measured winds from the mean wind values becomes insignificant „ 
Wind motions at these heights for this particular firing are 
characteristically nonturbulent „ 

Conclus ions 

The techniques of spectrum analysis hased on hydroc^ namic 
turbulence theory can be profitably applied to wind data obtained 
from sodium trail rocket firings, at least in the height range 
from 80 to 100 km„ The analysis of further firings should 
indicate whether or not the mean profiles determined as a prelude 
to spectrum analysis have any meteorological significance. 
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Appendix I 

Resume of the Computer Program 

The computer program has been written in FORTRAN IV 
for execution on an IBM 7094, and consists of the following 
routines : 

1, The main program SHEAR which reads and does some processing 
of the input data (which is punched on cards) , 

SHEAR calls subroutines 

a) FIT, which fits polynomial profiles to the zonal and 
meridional wind velocity/ height profiles by the method of least 
squares, utilizing 

b) MATS, to solve the set of independent linear equations 

c) The total wind speed profile is calculated by TRADE, 

d) OUTPUT performs the energy spectrum and correlation 
analyses, and printout of same, 

e) PAGE is a utility routine which turns and numbers the 
output pages, and prints an appropriate heading on each. 

Input Data 

The input data is on punched cards, and is read as follows: 
1, Header Card; the information punched on this card is 
printed at the top of each page of output. The full 72 columns 
can be used for any type of alphanumeric heading. 


2. The parameters specifying the mean zonal and meridional 
profiles o These determine the order of the polynomial fitted to 

a) the zonal (punched in columns 2,3) and 

b) the meridional (punched in columns 5,6) measured 
data to determine the mean profile . 

e.g., a 5 punched in column 3 will fit a polynomial of the form 

a^ + + a gh 2 + a 4^ 3 + a 5 ^ 

to the zonal wind data. 

The maximum allowable value of the subscript i in the a^ above 
is 10. 

3. The height range (must be integral, and in km) over which 
the energy spectrum and correlation analyses' are to be performed. 
This can be any fraction of the height interval covered by the input 
data, but must not be less than 10 km unless an appropriate 
spectral range is punched in column 13. The minumum of the height 
range interval is punched in columns 2 to 5, the maximum in 
columns 6 to 9. If OUTPUT is required to produce spectra over 
height differences to other than a maximum of 10 km, the required 
maximum height difference can be punched in columns 10 to 13. This 
entry must not exceed the height range specified in columns 2 to 9. 
If this last field is left blank, spectra to 10 km are output. 


4. The trial data, which must not exceed a maximum to 
minimum height range of 150 km„ The data is punched as year, 
month, day, hour (local time, 24 hour clock), height (km), 
irindspeed (metres/sec) , and azimuth (degrees) as follows 

YR MO DAY HOUR HEIGHT WINDSFEED AZIMUTH 


3 3 9 14 19 24 29 



all integer 

must be right adjusted if integer; 
anywhere in appropriate field if 
decimal punch is included. 

5. After all data of 4, a blank card to flag "end of data”. 

6. Further height ranges as for 3 as required. 

7. a) If a change in the mean wind profile is required a 
card with a 2 punched in column 1 may be substituted for any of 
the cards 6, and followed by a card of format 2 above which 
specified the new profile. This card is then followed by further 
height ranges (as for 6.) as required. 

7 b) If a completely new set of header plus data cards is 
to be processed after either stages 5, 6, or 7a, a card with a 1 
punched in column 1 will return control to the start of the program, 
which will then read in the new set of data cards sequenced as 


from 1 above 


APPENDIX II 


Listing of the FORTRAN IV Program 
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SHEAR 08/ 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT - INTERNAL 


SODIUM TRAIL HEIGHT SHEAR ANALYSIS PROGRAM 

CALCULATES THE ENERGY SPECTRUM FCR TOTAL SHEAR, AND FOR ZONAL ANC 
MERIDIONAL SHEAR COMPONENTS 

READS INPUT DATA AS FOLLOWS 

A HEACER CARD, PUNCHED WITH DETAILS OF FIRING TIME, ETC. 

FORMAT 12A6 

THE PARAMETERS SPECIFYING THE MEAN ZONAL AND MERIDIONAL PROFILES, 
FORMAT 213 

THE HEIGHT RANGE OVER WHICH ANALYSIS IS TO BE P ERF CRME C , ZM I N , 

ZMAX, ANC THE MAXIMUM OF THE OLTPUT SPECTRAL RANGE OESIRED. 

FORMAT 1X2F4.0, 14. 

SPECTRAL RANGE MUST NOT EXCEED ZMAX-ZMIN. 

IF NO SPECTRAL RANGE IS SPECIFIED, ZMAX-ZMIN MUST BE GREATER THAN 
10KM, ANC SPECTRA WILL eE OUTPUT TO IOKM. 

THE TRAIL DATA, YEAR MONTH CAY LOCAL TIME (HOURS AND MINUTES, 

24 HOUR CLOCK) HEIGHT(KM) W IND SPEE C ( METRE S / SEC ) WIND AZIMUTH 
( DEGREES ) . 

FORMAT 313, 1 5 , F 5 . 1 , 2F 5 . C 
A BLANK CARD 

FURTHER ZM IN, ZMAX AS REQUIRED. IF A 2 APPEARS IN COLUMN 1, 
PROGRAM WILL RE AC NEXT CARD AS NEW PROFILE SPECIFICATION, FOLLOWED 
BY FURTHER ZMIN, ZMAX CARDS. IF A 1 APPEARS IN COLUMN 1, THE 
PROGRAM WILL READ NEXT CARD AS HEADER CARD CF A COMPLETELY NEW SET 
OF DATA. 

DIMENSION eZ(lG), eM(lC) 

DIMENSION RESULT ( 1 2 ) , Z I ( 7 50 ) , W I ND I ( 750 ) , A ZR AD I ( 750 ) .ZONAL (750) , 
1EKI D ( 7 50 ) , TZCNAL(75C) ,T MERIC ( 75C ) 

DIMENSION TSPEEDI75C), AW(5) 

COMMON ZI, WINCI, AZRAC1, ZONAL, ERIC. TZONAL , TMERI D , T SPE ED 
C=0.05 
CS=U.0C00C1 
P I = 3 . 1415926 
TWOP 1 = 2 .0«P I 
999 NG0=0 
N I =0 
1=0 

NSUM=C 
ZER0=0 .0 

READ (2,1) RESULT 

1 FORMAT ( 12A6 ) 

17 READ (2,4)NP,NC 

4 FURMAT (213) 

NF I T = - 1 

2 REAU( 2 , 3 )NGO, ZMIN, ZMAX, NCIFF 

3 FORMAT ( I1.2F4.0, 14 ) 

IF ( NGC ) 19, 19, 18 

18 IF(NGC-2)999, 17, 17 

19 CONTINUE 
IF(NCIFF)20,20,21 

20 ND I FF= 10 

21 MIN=ZM IN+C 
MAX=7MAX+C 


SHEAR 08/ 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT - INTERNAL 


IF ( NI )5, 5,101 
TRAIL CATA INPUT 

5 READ l.?,6)MYEAR, MONTH, JOUR , LT IM E , Z , W I NO , AZDEG 

6 FORMAT ( 21 3 » I5»F5.1»2F5«C ) 

IF(MYEAR) 1CU» IOC ,7 

7 1=1+1 

I F ( 1-1) 11,11,8 

8 7DI FF=Z-Z It 1-1 ) 

M Z = 1 C . C «Z+C 

LZ= 1C .u »Z I ( I-1J+C 
LDI FF=M7-LZ 
I F ( LC I F F-2 ) 12,11,9 

12 WRITE (3,6)MYEAR»MCNTH»J0UR»LTIME,Z»WIND,AZDEG 
WR I TE ( 3 , l 3 ) 

13 FORMAT! 1X/////1X26FCATA CARCS OUT OF SECUENCE ////// 1X2CHE XECU T I CN 
1TERMINATEC//////1X ) 

PRINT 13 
CALL EXIT 

INTERPOLATION ROUTINE 

9 NPOINT=LC IFF/2-1 
NSUM=NSUM+NPOINT 

GRAO= ( W I NC - WINCI! I — 1 ) I/ZD1FF 
AZDIFF=AZCEG/57 . 3-AZRAC I ( 1-1 ) 

IF( ABSIAZCIFFI-P I 194,94,91 

91 IF! AZCIFF 192,94,93 

92 AZDIFF=AZCIFF+TWCPI 
G0T094 

93 AZDIFF=AZCIFF-TWCPI 

94 AZGRAC=AZCIFF/ZC IFF 
DU1QJ=1 ,NPO INT 

Z I ( I )=Z I ( 1-1 )+0.2 
W I NO I ( I ) = W I ND I ( 1-1 )+C.2*GRAC 
AZRACI ( I )=AZRAC I ( 1-1 )+C .2*AZGRAD 
I F ( AZRAL'I ! I ) 195,96,96 

95 AZRACI! I )=AZRACI (H+TWCPI 

96 IF! AZRACI ( I )-TWCFI )9e,9e,97 

97 AZRACI! I)=AZRACI! I )-TWCPI 

98 CONTINUE 
1=1 + 1 

10 CONTINUE 

11 ZI l I ) =Z 
WINCI ( I ) =W IND 
AZRACI! 1 )=AZDEC/57.3 
GO TC 5 

1G0 N I = I 

ALL DATA IN. COMMENCE CCMPLTATICN 
HD I F F = Z I(MI-ZKl) 

H P L U S = / I ( NI I+ZI ( 1 ) 

CALCULATE ZONAL ANC MERICICNAL WIND COMPONENTS 
DU2Q1 J= 1 , NI 

ZONAL 1 J )=WINDI ( J )• SIN! AZRACI ( J) ) 

201 ERID!J)=WINDI(J)« CCS l A ZR AC I ( J ) ) 

101 NZ=0 

lmin=iu*min 

LMAX=lu«MAX 


SI-EAR 08. 

EXTERNAL FCRMULA NUMBER - SOURCE STATEMENT - INTERNAL 

CALL P AGE ( RESULT , ZERO 
DO 10A J = 1 , N l 
LZ = Z I ( J )*1C.C + C 
IF(LZ~LMIN)10A» IC2» 1C2 

102 IF(LZ-LMAX)103» 1C 3 * 1C A 

103 N Z = N Z ♦ 1 
10A CONTINUE 

WRITE ( 3, 1C 5 )NI»NSUM»ZMIN,ZMAX ,NZ,NP,NC 

105 FORMAT ( 1X/////1X33HTCTAL NUMBER CF INPLT DATA POINTS I 8/// IX26H I NC 
1LUDING INTERPOL AT ICN OF 15, 7H POINTS / / / / 1 XA 5HNUMBER CF DATA POINT 
2S WITHIN THE H.EICH.T RANGE F5.C,6H KM TC,F5.C,3H KM I8////1X31HEAST 
3 WEST PROFILE SPECIFICATION 15 / / / 1 X3 3HN0RTH SOUTH PROFILE SPECIFI 
ACATICN 13 ///IX) 

S T ART = 1 

CALL PACE(RESULT, START ) 

WRITE (3,106 )ZMIN,ZMAX 

106 FORMAT ( 1X50HCALCLLATEC ZONAL AND MERIDIONAL MEAN WIND PROFILES 
110X12PHEIGFT RANGE F5.C,6H KM TO F5.C,2H KM/1H0 

2/10X 6HHEIGHT 1CX5HZCNAL 13X10HMERI CIONAL 1CX10HWIND SPEED /IX/ 

31X, 18X, 3IAX12HCATA MEAN, AX), /IX/) 

NHITE=(ZMAX-ZMIN)/2C.C-C 
NH I T E = NH I T E + 1 
I NK=5*NH I T F — 1 
N F I T = N F IT + 1 
IFINFIT 1202,202,203 

C CALCULATE COEFFICIENTS CF MEAN PROFILE POLYNOMIALS 

2C2 CALL NCRMAL I ZCNAL, ERIC,ZI ,NI,BZ,EM, NP,NG) 

203 CONTINUE 

C DETERMINE AND PRINT CUT MEAN WINC PROFILES 

SUDS = ( ZMAX+ZMIN-Z I(1))«5.G+1.0+C 
DU108L=MIN,MAX,NF ITE 
IZ=SUES-FLCAT I L *5 ) 

J=MAX+M IN-L 
H= J 

S=(2.0»H-HPLUS)/FC IFF+CS 
AZON A L = 0 • 0 
AMERI U=0.0 
DU20AM= 1 , 1C 

AZONAL-AZCNAL+DZ (M )*S*« (M-l ) 

AMERIC=AMERID+BM(M)*S*«(M-1) 

20 A CONTINUE 

MZONAL= AZCNAL 
MMERIC = AMER IC 

MSPEEC=SCRT ( AZCNAL«*2+AMERIC**2 ) 

WRI TE ( 3 , 107 )J, ZCNAL ( IZ ),MZCNAL*ERID I I Z ) , MMCR 1 0 , W I ND I ( I Z) , MS PEED 

107 FORMAT! 1X113, AX, 3( F 8 . C , AX I A , AX ) , / IX ) 

108 CONTINUE 

CALL PAGE ( RFSULT, START ) 

WRI TE ( 3 , 1C 9 )ZM IN ,ZMAX 

109 FORMAT ( IX , 22HTLRBULEN T WINC PRC F I LE , 10X 12HHE I GH T RANGE F5.0, 

1 t>HKM TC F5.0.3H KM / 1 X / IX , 1CX6H HE IG HT , 1C X6H T ZON AL , 1G X 1 1HTME R I D l CN A 
2L,10X6HTSPEED/1X) 

C DETERMINE AND PRINT CUT TURBULENT WINC PROFILES 

NH= ( MAX-M IN ) *5+1 
K=(ZMAX-ZI ( 1 ) )*5.C + 2.C + C 
NK= INK 


SHEAR 08 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT - INTEPNAL 

R M S Z = C .0 
RMSM=U.O 
RMS T = 0.0 
SUMSCT = 0 .C 
D01 12 J-= 1 * NH 
T ZONAL ( J ) =C .0 
TMER IC ( J ) =0 .0 
K = K-1 

H=ZMAX-U.2*FL0AT ( J — 1 ) 

I M = H + C 

S = ( 2.0*H-HPLUS ) /FCIFF+CS 
N K = N K + 1 
DU205M = 1 » 10 

TZONAH J)=TZONAL(J )+BZ(M)*S**(N-l) 

TMERIC! J) =T MERIC! J )+eN(M)*S«*(M-l) 

2C5 CONTINUE 

T ZONAL! J )= ZONAL (K)-TZCNAL(J) 

TMERIC!J)=ERIC(K )-TMERIC< J ) 

CALL TRACE ( J » K ) 

C CALCULATE RMS TURBULFNT VELOCITIES FOR THE HEIGHT RANGE ZMIN/ZMAX 

RMSZ=TZCNAL ( J ) * « 2 + RM S Z 
RMSM = T MERIC ( J )«*2+RMSM 
RMST=TSPEEC I J )**2+RMST 
SUMSCT=SUMSGT+1.C 
I F { NK- INK ) 112, 1 12, 1 1C 

110 WRITC! 3, 111 ) IH, TZCNAL! J ),TMERIC( J ), TSPEED ( J ) 

111 FORMAT I 1X» 10X14, 12XFf.C»llXF8.C* 13XF5.C/1X) 

NK = 0 

112 CONTINUE 

CALL PAGEIRESULT, START) 

RMSZ = SCRT (RMSZ/SLMSCT ) 

RMSM= SC RTIRMSM/ SUMSCT ) 

RMS T = SCRT ( RMS T/ SUMSCT ) 

WRITE! 3,206 ) ZM IN, ZM AX, < BZ ( I ), 1 = 1, NP ) 

206 FORMAT! 1X1CHNCRMAL IZEC 

1 1X12 HH EIC-HT RANGE10XF7.C, 7H KM. TO F7.0,4H KM./1H0/ 

2 1X42HC0EFFIC 1ENTS GIVING BEST FIT TO ZONAL C ATA /I X/ 10E 1 1 . 3 ) 
WRI TE(3,207 ) (BM( I ) , I=1,NC ) 

207 FORMAT! IX/ IX/ 

11X/ 1X47HCCEFFIC IENTS GIVING BEST FIT TC MERIDIONAL D A T A / 1 X/ 1 X ICE 1 1 
2.3/1X/) 

WRI TE( 3,208 )RMSZ ,RMSM,RMST 

208 FORMAT ( IX/ IX 

l 1X22HRMS TURBULENT VELGC I T Y / 1 X/ I X 1C X5H ZON ALF 1 1 . 0 , 1 1 HME TRE S/ S 

2EC. / 1X/1XIGX10HMER IC IONAL F 6 . C , 1 1 HM ETR E S / SEC . / 1 X / 1 X 10X 10HWI NO SPEE 
30 F6.C, 11HMETRES /SEC. ) 

CALL PAGEIRESULT, START ) 

WRITE (3,113)ZMJN»ZMAX 

113 FORMAT! lXAtHENERGY SPECTRUM OF ZCNAL WIND HEIGHT VARIATION, 

110 X 12 H HEIGHT RANGE FS.C,6H KM TO F5.0,3H KM/1X) 

C CALCULATE AND OUTPUT ENERGY SPECTRUM FUNCTIONS 

CALL CLTPUT ! TZCNAL ,NH,NCIFF ) 

CALL PAGE(RESULT, START ) 

WRITE (3,114) ZM1N.ZMAX 

114 FORMAT! IX51HENERGY SPECTRUM OF MERIDIONAL WIND HEIGHT VARIATION, 



o o o 


/ $ 1 0 J02T RGR ELCC- 11 

(PAUSE 

(EXECUTE IEJ08 

IBJCB VERSION 2, 7090 - PR - 929 

1 


0 

(IB JOB GC 

(IBFTC OUT M9A,XR7 

1 CUT 

EXTERNAL FORMULA NUMBER 


0 


SUBROUTINE OUTPUT ( T SPEEC , NH,ND IFF ) 


SOURCE STATEMENT 


INTER! 


PRODUCES HEIGHT SPECTRUM AND PERFORMS CORRELATION ANALYSIS. 


DIMENSION TSPEEC (75C ) 

WRITE ( 3 • 1 ) 

1 FORMAT! IX, 10X7HCELTA H , 5X 10FE I CELTA H), 
15X11HCCRRELATICN»3X13HS IGNIF 1-G /IX) 

NEND = 5 «NC IFF 
L I NE=0 

D0125K* 1 , NEND 

L I NE*L I NE+ 1 

DELTAH=FLCAT(K)«C.2 

SUM = 0 .0 

ENERGY*0 .0 

ERGX*0 .0 

GJPKX J=0.0 

SQ J*0 .0 

SQ JPK*0 .0 

D0123 J= 1,NH 

IF(NH-J-K) 124, 122, 122 

122 JPK= J+K 

ENERGY* ENERGY*! TSPE EC! JPK J-TSPEED1J ))*«2 
SUM=SUM+1.0 

GJPKX J*GJPKXJ+TSPEECUFK)*T SPEED! J) 
SQJ=SCJ+TSPEED(J)**2 
SQJPK*SQJPK+T SPEED! JFK)««2 

123 CONTINUE 

12A IF1SUM ) 117, 117, 118 

117 G = 0 .0 

S IGNIF*9.99 
G0T0119 

118 ENERGY* ENERGY /SUM 
G=GJPKXJ/SCRT!SQJ*SCJPK ) 

S IGNIF* <1.0-G*«2) /SORT! SUM-1.0) 

119 DI FF= 1 .0-G 
IF1LINE-5C )125, 125, 116 

116 WRITE ! 3 ,2 ) 

2 FORMAT! 1H1/1X) 

WRITE (3*1) 

L I NE=Q 

125 WRITE! 3 ,115 JOELTAH, ENERGY, G,SI GN IF,D IFF 

115 FORMAT! IX, 6XFB. 1 , 7XF8.1 , 8XF8. 3 ,F 11 . 3 ,F8. 2 ) 

RETURN 
END 



SHEAR 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT 

UOX12HHLIGHT RANGE F5.0.6H KM TC F5.0.3H KM/1X) 

CALL CUTPLTITMERID.NF.NCIFF) 

CALL PAGE( RESULT , START ) 

WRITE ( 3, 115 IZMIN.7MAX 

115 FORMAT ( 1X46HENERGY SPECTRUM OF WIND SPEED HEIGHT VARIATION 
110X12H HEIGHT RANGE F5.C.6H KM TO F5.0.3H KM/1X) 

CALL OUTPUT (TSPEEC»NH,NCIFF ) 

GOT02 
EN 


08 

INTERNAL 


1 


$ I B JOB GC 

$ l BF TC NORMAL M94,XR7 


NORMAL 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT - INTE 


SUBROUTINE NORMAL ( ZC N A l , ER ID , 1 1 , N I , B Z , EM , N P , NQ ) 

DETERMINES BEST FIT PROFILE TO THE NORMAL 1 7EC TOTAL HEIGFiT RANG 

DIMENSION ZONAL! 75C ) , ER 10(750), ZII750) * B Z ( 1C) »BM( 1C) , ZIN0RM(750 
D01 J = 1 1 N I 

Z I NORM! J)*(2.0*Z I! J )-ZI(M)-ZI(l))/(ZI(M)-ZI(l)> 

1 CONTINUE 

CALL FIT (ZONAL. NP , BZ , N I , Z INORM ) 

CALL FIT (ERIC, NO , EM , N I , Z I NORM ) 

RETURN 



03 / 


$IBJCB GC 

$IBFTC TRACE f'S's.Xrt? 
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$ I BFTC FIT 


M V 4 1 X R 7 


o o o 


FIT 08/17 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT - INTERNAL FO 


SUBROUTINE FIT I BLOW , N S PEC , A , N , X ) 

FITS A HEIGHT PROFILE SPECIFIED BY NSPEC TO THE MEASURED DATA 
BLOW, DETERMINING THE COEFFICIENTS A BY THE METHOD OF LEAST 
SUUARES. PROFILE FITTED OVER TOTAL HEIGHT RANGE OF INPUT DATA. 

DIMENSION BLOW (750), A(10), X(7 SO), S(1C,11) 

0031=1,10 

3 A l J ) = 0 • 0 

IF ( NSPEC )l :.'0, 100, A 

4 NADD=NSPEC+1 
DUlOJ=l, NSPEC 
DU10K= 1 , N A D D 

10 S ( J , K ) =0 • U 
1)05 J= 1 , NSPEC 
D05K=1 , 4SPEC 
DJ5I=1»N 

5 S(J,K)=S(J»K )+X( I ) * * ( K + J- 2 ) 

DUfaJ=l, NSPEC 

K = 0 

D U 6 I = 1 , N 
K = K+1 

6 S(J,NADD)=S(J,NACD)+X( I )*»(J-1)«6L0W(K) 

CALL MATS! S,A,NSPEC, V ISS) 

IF ( MISS) 100, 100,7 

7 WRITE! 3,8) 

8 FOP. MAT ( 1X75 HERR OR IN INPUT DATA HAS RESULTED IN MATRIX S BEING UNS 
1UITABLE FOR INVERSION / / // 1X20HEXEC UT I C.N TERMINATED ///////) 

PRINT 3 
CALL EXIT 
100 RETURN 


o n o 


MATS 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT 


SUBROUTINE MATS ( S , A, NSPEC , M I S S ) 

REDUCES THE AUGMENTED MATRIX S TO PRODUCE THE COEFFICIENTS 

DIMENSION S( ID, 1 1 ) . A ( 1 0 ) 

M l S 5 = - l 
MM=NSPEC+1 
N=NSPEC 
DU 1 5 I = 2 f N 
70 11=1-1 

7 DO 1 5 J= 1 1 I I 

8 I F ( S ( I , J ) >9,15,0 

9 IF ( AESISI J, J ) ) - APS ( SC I, J ) ) ) 11, 10,10 

10 R=S( I , J ) / S ( J » J ) 

GOTO 1 3u 

11 P.= S( J, J)/S( I, J> 

D012K=1,MM 

B = S( J,i<) 

S(J»K)=S(I,Kj 

12 S ( I , K ) = B 
130 JJ=J+1 

13 DJ14K= J J ,MM 

14 S ( I » K ) =S ( I » K ) -R*S ( J » K ) 

15 CONTINUE 

IF ( ARS( S(N,N) )-l .OE— 10 ) 16, 16, 17 

16 MI SS = + 1 

G0T029 

17 A(N)=S(N,MM)/S(N,N) 

0028 1=2, N 
JJ=N-I+1 

B = 0 . 0 
I I = N— I +2 
D 0 2 5 K = I I , N 
2b B=B+S( JJ,K)*A(K) 

IF(ABS(S(JJ,JJ))-1.0E-10)16,16,28 

28 A(JJ)=(S(JJ*MM)-B)/S(JJ»JJ) 

29 RETURN 
END 


08/17 
INTERNAL FO 


A 
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ilBFTC PAGE M94,X*7 


o o o n 


PAGE 

EXTERNAL FORMULA NUMBER - SOURCE STATEMENT 


SUBROUTINE PAGE ( RESULT, START ) 

TURNS PAGE, NUMBERS IT, AND WRITES HEADING AS APPEARING 
ON RESULT CARD. 

DIMENSION RESULT (12) 

IF(START) 1,1,2 

1 N0P=0 

2 NOP=NOP+l 

WRITE { 3, 3) RESULT, NOP 

3 FORMAT! 1H1 / IX 12 A 6 , 30X AH PAGE 1 5/ / ) 

RETURN 


08/17 
INTERNAL FO 


APPENDIX III 


SAMPLE DATA AND COMPUTER OUTPUT 


Noc ! 



63 

5" 

21 

1910105.6 

140" 

130 

6 3 

5 

21 

1910103.8 

130 

1*28 

63 

5 

21 

1910107.0 

125 

132 

63 

3 

21 

1910108.0 

92 

140 

63 

5 

21 

1910109.0 

50 

160 

6 3 

5 

21 

1910110.0 

42 

240 

63 

T" 

TT 

1910111.0 

56 

250 

63 

3 

21 

1910113.0 

60 

262 

63 

3 

21 

1910114.0 

84 

274 

63 

5 

21 

1910115.0 

68 

270 

63 

3 

21 

1910117.2 

72 

282 

63 

3 

Ti“ 

T910125'.0 

43 

286 


63 5 21 1910130.0 30 296 

63 5 21 1910133.0 14 310 

63 5 21 191 0140 75 15 316” 


— BLANK, CARl>~ 



PRELIMINARY PAPER EGLIN DATA* MAY 21* 
5" 5 * 

8'-; 1PL 

~63 WTI'910 "69.0 34 256 

63 " T 2TT9T6 70.6 31 256 

“63 - TTl lVlcTfiTo 2 7 270~ 

" 6 3~ 5 2 1 " 1 9 1C” 72 '.~4 4 3~ “ 212 

" '6 r "5'21 1910“T2V6 2 o" ~2 90 

~63 5 2T"l 9'uT 77~."o 41 256 

63 5 21 1910 80. 4 25“ 358' 


63 

— T 

21 

1910 82.0 

08 

310 

63 

5 

21 

1910 82.4 

28 

202 

63 

3 

21 

1910 83.6 

27 

200 

63 

3 

21 

1910 84.0 

54 

190 

63 

"3 

21 

i 9 1 0 ~83 .0 

’ 3 4 

194 

63 

3 

21 

1910 “3 6..~4~ 

30 

108 

63 

5 

21 

1910 88.0 

60 

0 80 

63 

3 

21 

1910 91.6 

76 

093 

63 

'5 

21' 

1910 93.0 

30 

178 

6^ 

— 5 - 

2T 

'1^1T7'96.(T 

40 

" 120 

6 3 

3 

21 

1910 96. 6 

54 

130 

63 

3 

21 

1910 97.6 

4 0 

115 

63 


TT 

1910 99 TO 

‘ '42 

■‘"084 

6 3 

3 

21 

1910 99.6 

" 64 

07S 

63 

5 

21 

1910100.0 " 

3 8 

070 


“”3 

yr 

1 91 ofOoTfl 

~~ 65 ~ 

069 

63 

5 

21 

1910102.0 

8 0 

079 

63 

3 

21 

1910103.0 " 

93 

094 

63 

3 

“21" 

iTT 0104TT) 

nr 

104 

63 

3 

21 

1910104.2 

in 

120 

63 

y 

21 

1910105.0 

12 9 

122 


I 


PROF I LE » 5 * 5 


1963. 16/9/64. 




TRIAL ANALYSIS. EGLIN DATA, M AY 21, 1963. 15/9/64 


PROFILE, 5, 5. 


TOTAL NUMBER OF INPUT CAT* POINTS 356 
INCLUCING INTERPOLATION OF 313 PGINTS 

NUMBER OF CATA POINTS WITHIN THE HEIGHT RANGE 80. KM TO 100. KM 
EAST VEST PROFILE SPECIFICATION 5 


101 


NORTH SOUTH PROFILE SPECIFICATION 5 


TRIAL ANALYSIS. EGLIN DATA, MAY 21, 1963. CONTINUOUS PROFILE, 5, 5. 
CALCULATED ZONAL AND MERIDIONAL MEAN WIND PROFILES HEIGHT RANGE 


HEIGHT 

ZONAL 


MERIDIONAL 

WIND 

SPEED 


DATA 

MEAN 

DATA 

MEAN 

DATA 

MEA' 

100 

TO. 

34 

20. 

-21 

72. 

40 

9s 

56 . 

37 

21 . 

-21 

60. 

4 3 

98 

42. 

39 

4. 

-22 

42. 

45 

9 7 

39. 

41 

-11. 

-22 

41. 

47 

96 

34. 

43 

-35. 

-22 

48. 

49 

9 b 

31 . 

45 

-23. 

-22 

42. 

50 

94 

6 . 

46 

-51. 

-22 

52. 

51 

9 3 

31 . 

46 

-51 . 

-21 

59. 

51 

9 2 

5 5. 

47 

-38. 

-21 

67. 

51 

9 1 

73. 

47 

-13. 

-21 

74. 

51 

90 

72. 

46 

-1. 

-2 J 

72. 

50 

39 

67. 

45 

4. 

-19 

67. 

49 

88 

62. 

43 

8. 

-18 

63. 

47 

8 7 

55. 

41 

5. 

-18 

55. 

45 

86 

44 . 

39 

-4. 

-16 

44. 

43 

3 u 

31 . 

36 

-9. 

-15 

32. 

39 

34 

-12. 

33 

-53. 

-14 

54. 

36 

8 3 

-9. 

29 

-25. 

-13 

27. 

32 

32 

-10. 

24 

-26. 

-11 

28. 

27 

31 

-8. 

19 

10. 

-10 

12. 

22 

8 0 

-3. 

14 

23. 

-8 

23. 

Id 


TRIAL ANALYSIS. EGLIN DATA, MAY 21, 1963 


CONTINUOUS PROFILE, 5, 5. 


TURBULENT WIND PROFILE HEIGHT RANGE 30. KM TO 100. KM 


HEIGHT 

TZONAL 

TMERIDIONAL 

TSPEED 

100 

35. 

41. 

42. 

99 

19. 

42. 

26. 

99 

2. 

26. 

16. 

97 

-3. 

11 . 

31. 

96 

-1C, 

-12. 

32. 

95 

-14. 

-6. 

29. 

94 

-4 0. 

-29. 

3. 

9 3 

-16. 

-29. 

31. 

92 

3. 

-16. 

53. 

91 

26. 

8. 

62. 

9 0 

25. 

19. 

52. 

S9 

22. 

24. 

45. 

88 

1 8 . 

27. 

38. 

87 

13. 

24. 

33. 

86 

4. 

13. 

30. 

8 5 

-6. 

7. 

27. 

84 

-4 5. 

-38. 

18. 

53 

-3 9, 

-12. 

-5. 

82 

-35. 

-14. 

0. 

81 

-28. 

20. 

3. 

3 0 

-18. 

31 . 

10. 



TRIAL ANALYSIS. ECLIN CATA, MAY 21, 1962. 15/9/64. 

NORMALIZED HEIGHT RANGE 8C. KM. TC ICO. KM 

COEFFICIENTS GIVING BEST FIT TC ZCNAL CATA 
0.386E 02 -C.126E 03 -0.244E C3 0.174E C3 C.205E 03 

COEFFICIENTS GIVING GEST FIT TC MERICICNAL CATA 
-0.304E 02 C.198F 02 C.148E 03 -0.168E 02 -0.125E 03 

RMS TURBULENT VELOCITY 

ZCNAL 22 .METRES/SEC. 

MERICICNAL 22 .METRES/SEC. 


PRCF I LE , 5 , 5 . 


KIND SPEED 17 .METRES/SEC 


TRIAL ANALYSIS. EGLIN DATA, MAY 21, 1963 


CONTINUOUS PROFILE, 5, 5. 


ENERGY SPECTRUM OF ZONAL WIND HEIGHT VARIATION HEIGHT RANGE 80. 


DELTA H 

DELTA V ** 2 

CORRELATION 

SI GNIF 

1-G 

0.2 

24.3 

0.979 

0.004 

0.02 

0.4 

64.9 

0.943 

0.011 

0.06 

0.6 

113.7 

0.900 

0. 019 

0.10 

0.8 

169.4 

0.851 

0. 028 

0.15 

1.0 

230.7 

0.797 

0.037 

0.20 

1.2 

297.1 

0.739 

0. 047 

0.26 

1.4 

369.4 

0.676 

0. 056 

0.32 

1.6 

450.4 

0 .604 

0. 066 

0.40 

1 .8 

534.8 

0.525 

C. 076 

0.47 

2.0 

619.9 

0 .449 

0. 084 

0.55 

2.2 

70 3.2 

0.374 

0. 091 

0.63 

2.4 

78 9.0 

0 .297 

0. 097 

0.70 

2.6 

877.1 

0.217 

0. 102 

0.78 

2.e 

964.9 

0.136 

0. 106 

0.86 

3.0 

1054.0 

0 .053 

0. 108 

0.95 

3.2 

1133.4 

-0.023 

0. 109 

1.02 

3.4 

1209.9 

-0.099 

0. 109 

1.10 

3.6 

1285.2 

-0 .176 

0. 107 

1.18 

3.8 

1360.8 

-0.257 

C. 194 

1.26 

4 . 0 

1436.0 

-0.342 

0.099 

1.34 

4.2 

1507.0 

-0.429 

0.092 

1.43 

4.4 

1567.1 

-0 .514 

0. 083 

1.51 

4.6 

1614.3 

-0 .597 

0. 073 

1.60 

4.8 

1681.4 

-0.651 

0. 066 

1.65 

5.0 

1745.7 

-0 .696 

0. 060 

1.7C 

5.2 

1304.0 

-0.735 

0. 053 

1.74 

5.4 

1852.3 

-0.766 

0. 048 

1.77 

5.6 

1885.5 

-0 .788 

0. 045 

1.79 

5 • 8 

1897.3 

-0.900 

0. 043 

1.80 

o . 0 

1980.3 

-0 .301 

0.043 

1.30 

6.2 

1863.6 

-0 .795 

0. 044 

1.80 

6.4 

1847. 4 

-0 .783 

0. 047 

1.78 

6.6 

1827.7 

-0.762 

0. 051 

1.76 

6.8 

1803.3 

-0.730 

0.057 

1.73 

7.0 

1775.0 

-0.690 

0.065 

1.69 

7.2 

1740.9 

- 0.642 

0. 073 

1.64 

7.4 

1701.2 

-0.588 

0. 082 

1.59 

7.6 

1656.5 

-0.528 

0. 092 

1.53 

7.8 

1608.4 

-0 .466 

C. 100 

1.47 

8.0 

1555.6 

-0.401 

0. 198 

1.40 

8.2 

1497.1 

-9.333 

0. 116 

1.33 

8.4 

1432.6 

-0.263 

0. 122 

1.26 

8.6 

1361.9 

-0.190 

0. 128 

1 .19 

8.8 

1285.1 

-0.116 

0. 132 

1.12 

9.0 

1202.6 

-9.040 

0. 135 

1.04 

9.2 

1114.7 

0.939 

G. 136 

0.96 

9.4 

1020.5 

0.120 

0. 135 

0.8 8 

9.6 

920.9 

0.207 

0. 133 

0.79 

9.8 

819.6 

0 .296 

C. 128 

0.70 

10.0 

733.6 

0 .372 

0. 122 

9.63 



TRIAL ANALYSIS. EGLIN DATA, MAY 21, 1963. CONTINUOUS PROFILE, 5, 5. 
ENERGY SPECTRUM OF MERIDIONAL KINO HEIGHT VARIATION HEIGHT RANGE 


DELTA H 

DELTA V»*2 

CORRELATION 

SI GNIF 

1-G 

0.2 

31.3 

0.973 

0.005 

0.03 

0.4 

94.1 

0.917 

0.016 

0.08 

0.6 

164.4 

0.853 

0.028 

0.15 

0 • 8 

241.7 

0.781 

0.040 

0.22 

1.0 

326.3 

0.701 

0.052 

0.30 

1.2 

418.3 

0.612 

0. 064 

0.39 

1.4 

615.2 

0.518 

0. 076 

0.48 

1.6 

616.0 

0 .420 

0.036 

0.58 

1.8 

723.6 

0.316 

0. 094 

0.68 

2.0 

825.3 

0.219 

0. 100 

0.78 

2.2 

935.1 

0.116 

0. 105 

0.88 

2.4 

1C52.2 

0 .00 8 

0. 107 

0.99 

2.6 

1175.9 

-0.104 

0. 106 

1.10 

2.8 

1305.7 

-0.219 

0. 103 

1.22 

3 . 0 

1429.3 

-0.325 

0. 097 

1.33 

3.2 

1534.6 

-0.411 

0.091 

1.41 

3.4 

1614.7 

-0 .479 

0. 085 

1.48 

3.6 

1667.7 

-0.535 

0.079 

1.54 

3.6 

1711.0 

-0.582 

0.074 

1.58 

<t.C 

1 741 .1 

-0 .617 

0.069 

1.62 

4.2 

1756.3 

-0 .641 

0.066 

1.64 

4.4 

1758.3 

-0.656 

0. 065 

1.66 

4.6 

1745.6 

-0 .662 

0. 064 

1.66 

4.8 

1776.6 

-0 .680 

0. 062 

1.68 

3.0 

1836.6 

-0.715 

0. 056 

1.72 

5.2 

1878.5 

-0.733 

0.054 

1.73 

3.4 

190U.7 

- 0.735 

0. 054 

1.74 

5.6 

1902.8 

-0 .723 

0.056 

1.72 

5.8 

1885.6 

-0 .699 

0.061 

1.70 

6 • 0 

1853.4 

-0 .666 

0. 066 

1.67 

6.2 

1809.8 

-0 .625 

0. 073 

1.63 

6.4 

1751.9 

-0.573 

0.031 

1.57 

6.6 

1679.1 

-0 .510 

0. 090 

1.51 

6.8 

1592.6 

-0.435 

0. 100 

1.44 

7.0 

1494.3 

-0.350 

0. 109 

1.35 

7.2 

1385.9 

-0 .257 

0. 117 

1.26 

7.4 

1269.8 

-0 .156 

0. 123 

1.16 

7.6 

1 149.0 

-0 .051 

0. 127 

1.05 

7.8 

1028.4 

0.057 

0. 128 

0.94 

3 . C 

912.5 

0.164 

0. 126 

0.84 

3.2 

802.4 

0.266 

0. 121 

0.73 

3.4 

697.1 

0 .366 

0. 114 

0.63 

8.6 

598.6 

0 .459 

0. 105 

0.54 

8.8 

50e.4 

0.544 

0.0 94 

0.46 

1 . 0 

427.7 

0 .620 

0. 083 

0.38 

9.2 

357.9 

0.684 

0. 092 

0.32 

9.4 

302.8 

0. 735 

0. 063 

0.27 

9.6 

260.1 

0 .774 

0.056 

0.23 

9.3 

231.1 

0.301 

0. 050 

0.20 

10.0 

219.4 

0.812 

0. 048 

0.19 



TRIAL ANALYSIS. RGLIN DATA, MAY 21, 1963. CONTINUOUS PROFILE, 5, 5. 


ENERGY SPECTRUM OF WIND SPEED HEIGHT VARIATION HEIGHT RANGE 


DELTA H 

DELTA V**2 

CORRELATION 

SIGNIF 

1-G 

0.2 

15.1 

0.993 

0.001 

0.01 

0.4 

47.1 

0.979 

0.004 

0.02 

0.6 

83.2 

0.964 

0.007 

0.04 

0.8 

122.0 

0.947 

0.011 

0.05 

1.0 

161.5 

0.930 

0.014 

9.07 

1.2 

200.6 

0.914 

0.017 

0.09 

1.4 

240.7 

0.898 

0.020 

0.10 

1.6 

286.9 

0.879 

0.024 

0.12 

1.8 

337.0 

0.859 

0.028 

0.14 

2.0 

382.2 

0.841 

0.031 

0.16 

2.2 

421.0 

0.827 

0.034 

0.17 

2.4 

455.3 

0.815 

0.036 

0.19 

2.6 

483.4 

0.805 

0.038 

0.20 

2.8 

502.9 

0 .799 

0.039 

0.20 

3.0 

512.5 

0.797 

0.040 

0.20 

3.2 

517.3 

0.796 

0.040 

0.20 

3.4 

523.9 

0.793 

0.041 

0.21 

3.6 

542.8 

0 .789 

0.042 

0.21 

3.3 

561.7 

0.783 

0. 043 

0.22 

4.0 

584.0 

0.775 

0.045 

0.22 

4.2 

60 6.6 

0.768 

0.046 

0.23 

4.4 

628.3 

0.761 

0.048 

0.24 

4.6 

650.9 

0.755 

0.049 

0.25 

4.8 

679.5 

0.746 

C. 051 

0.25 

5.0 

712.5 

0.734 

0.053 

0.27 

5.2 

748.8 

0.722 

0.056 

0.28 

5.4 

785.1 

0.710 

0.058 

0.29 

5.6 

818.7 

0.700 

0. 060 

0.30 

5.8 

847.5 

0.692 

0.062 

0.31 

6.0 

370.6 

0.686 

0.063 

0.31 

6.2 

894.8 

0.680 

0. 065 

0.32 

6.4 

919.4 

0.674 

0. 066 

0.33 

6.6 

943.4 

0 .668 

0. 068 

0.33 

6.8 

968.1 

0.661 

0.069 

0.34 

7.0 

996.2 

0.653 

0.071 

0.35 

7.2 

1023.3 

0 .645 

0.073 

0.35 

7.4 

1048.3 

0.638 

0. 075 

0.36 

7.6 

1070.2 

0.630 

0. 077 

0.37 

7.3 

1083. 9 

0.623 

0. 078 

0.38 

8.0 

1103.4 

0 .616 

0.080 

0.38 

8.2 

1114.0 

0.608 

0. 082 

0.39 

8.4 

1121.1 

0.600 

0.084 

0.40 

8 . 6 

1120.2 

0.593 

0.036 

0.41 

8.8 

1114.6 

0.588 

0.087 

0.41 

9.0 

1103.8 

0.585 

0.089 

0.42 

9.2 

1C87.7 

0.584 

0.090 

0.42 

9.4 

1066.1 

0.586 

0,090 

0.41 

9.6 

1039.0 

0.590 

0. 090 

0.41 

9.8 

1007.8 

0.596 

0. 090 

0.40 

10.0 

973.4 

0 .603 

0.090 

0.40 


